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Abstract 

We derive a local-time path-integral representation for a generic one-dimensional time-independent 
system. In particular, we show how to rephrase the matrix elements of the Bloch density matrix as 
a path integral over ^-dependent local-time proHles. The latter quantify the time that the sample 
paths x(t) in the Feynman path integral spend in the vicinity of an arbitrary point x. Generalization 
of the local-time representation that includes arbitrary functionals of the local time is also provided. 
We argue that the results obtained represent a powerful alternative to the traditional Feynman-Kac 
formula, particularly in the high and low temperature regimes. To illustrate this point, we apply 
our local-time representation to analyze the asymptotic behavior of the Bloch density matrix at low 
temperatures. Further salient issues, such as connections with the Sturm-Liouville theory and the 
Rayleigh-Ritz variational principle are also discussed. 


I. INTRODUCTION 

The path integral (PI) has been used in quantum physics since the revolutionary work of Feynman [T] , although the 
basic observation goes back to Dirac Hi who appreciated the role of the Lagrangian in short-time evolution of the 
wave function, and even suggested the time-slicing procedure for finite, i.e., non-infinitesimal, time lags. Since then 
the PI approach yielded invaluable insights into the structure of quantum theory [3] and provided a viable alternative 
to the traditional operator-formalism-based canonical quantization. During the second half of the 20th century, the 
PI became a standard tool in quantum held theory [S] and statistical physics [5], often providing the easiest route 
to derivation of perturbative expansions and serving as an excellent framework for (both numerical and analytical) 
non-perturbative analysis [?]■ 

Feynman PI has its counterpart in pure mathematics, namely, in the theory of continuous-time stochastic pro¬ 
cesses [5] . There the concept of integration over a space of continuous functions (so-called huctuating paths or sample 
paths) had been introduced by Wiener [9] already in 1920’s in order to represent and quantify the Brownian motion. 
Interestingly enough, this so-called Wiener integral (or integral with respect to Wiener measure) was formulated 2 
years before the discovery of the Schrddinger equation and 25 years before Feynman’s PI formulation. 

The local time for a Brownian particle (in some literature also called sojourn time) has been of interest to physicists 
and mathematicians, since the seminal work of Paul Levy in 1930’s m- In its essence, the local time characterizes the 
time that a sample trajectory x{t) of a given stochastic process spends in the vicinity of an arbitrary point X. This 
in turn dehnes a sample trajectory of a new stochastic process. A rich theory has been developed for local-time 
processes that stem from diffusion processes (see, e.g., Ref. [12] and citations therein). For later convenience, we 
should particularly highlight the Ray-Knight theorem which states that the local time of the Wiener process can be 
expressed in terms of the squared Bessel process HMg. In contrast to mathematics, the concept of the local time is 
not uniquely settled in physics literature. Various authors define essentially the same quantity under different names 
(local time, occupation time, traversal time, etc.), and with different applications in mind. For example, in Ref. |16j 
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the traversal time is used to study quantum scattering and tunneling processes, in m the small-temperature behavior 
of the equilibrium density matrix is analyzed with a help of the occupation time, while in Ref. |18j the large-time 
behavior of path integrals that contain functionals of the local time is discussed. 

The aim of this paper is to derive a local-time PI representation of the Bloch density matrix, i.e., the matrix 
elements {xb\ e~^^ ja^a) of the Gibbs operator. This can serve not only as a viable alternative to the commonly used 
Feynman-Kac representation but also as a powerful tool for extracting both large and small-temperature behavior. 
Apart from the general theoretical outline, our primary focus here will be on the low-temperature behavior which is 
technically more challenging than the large-temperature regime. In fact, the large-temperature expansion was already 
treated in some detail in our previous paper m- Last but not least, we also wish to promote the concept of the local 
time which is not yet sufficiently well known among the path-integral practitioners. 

The structure of the paper is as follows. To set the stage we recall in the next section some fundamentals from the 
Feynman PI which will be needed in later sections. In Section [ITT| we provide motivation for the introduction of a local 
time, and construct a heuristic version of the local-time representation of PPs. The key technical part of the article 
is contained in Section |IV[ where we derive by means of the replica trick the local-time representation of the Bloch 
density matrix. Relation to the Sturm-Liouville theory is also highlighted in this context. A local-time analog of the 
Feynman-Matthews-Salam formula muD] is presented in Section]^ and its usage is illustrated with a computation of 
the one-point distribution of the local time. Since a natural arena for local-time Pi’s is in thermally extremal regimes, 
we confine our attention in Section |VI| on large- and small-/3 asymptotic behavior of the Bloch density matrix. 
There we also derive an explicit leading-order behavior in large-/3 (i.e., low-temperature) expansion. The analysis 
is substantially streamlined by using the Rayleigh-Ritz variational principle. Finally, Section m summarizes our 
results and discusses possible extensions, applications, and future developments of the present work. For the reader’s 
convenience the paper is supplemented with appendix which clarifies some finer technical details. 


II. PATH-INTEGRAL REPRESENTATION OF THE BLOCH DENSITY MATRIX 


Consider a non-relativistic one-dimensional quantum-mechanical system described by a time-independent Hamilto- 
V(x) where p |a;) = —ihdx |a:). Throughout this paper we will study the matrix elements 


nian H = 


p{Xa,Xb,P) = {Xb\e \Xa) 


( 1 ) 


of the Gibbs operator e , where /3 = l/{kBT) is the inverse temperature and feg is the Boltzmann constant. The 
matrix p(Xa,Xb, /3), known also as the Bloch density matrix, is a fundamental object in quantum statistical physics, 
as the expectation value of an operator O at the temperature T can be written in the form 


( 6 ) = 


1 


dXadXb p{Xa, Xb, jd) {xb\0\Xa) , 


( 2 ) 


where Z = J^dx p{x, x, fd) is the partition function of the system. In case of need, ensuing quantum mechanical 
transition amplitudes can be obtained from Q via a Wick rotation which formally amounts to the substitution 
/3 —>■ it/h, converting thus the Gibbs operator to the quantum evolution operator 

The matrix elements in Eq.Q can be represented via the path integral as HIT] 


p{Xa,Xb,P) 


x{^h)—xb 


x{0)—Xa 


T>x{t) exp 



'M 

Y 




V{x) 


(3) 


This represents a “sum” over all continuous trajectories x{t), t G [0,fdh], connecting the initial point a;(0) = Xa with 
the final point x{(dh) = Xb- It should be noted that the integral J^^dr + V{x)\ is the classical Euclidean action 
integral along the path x{t) with 0 < r < [dh. In the following we will denote the Euclidean action as A. The integrand 
in A, i.e., {M/2)x^{t) -|- V{x{t)), can be identified with the classical Hamiltonian function, in which the momentum 
p is substituted for Mx. One can also regard (bl as an expectation value of the functional exp[— dr V{x{T))/h] 
over the (driftless) Brownian motion with the diffusion coefficient M/2h, and duration fdh, that starts at point Xa, 
and terminates at Xb- The latter stochastic process is also known as Brownian bridge. 


HI. LOCAL-TIME REPRESENTATION OF PATH INTEGRALS: HEURISTIC APPROACH 

The purpose of this section is twofold. Firstly, we would like to motivate a need for reformulation of Pi’s in the 
language of local-time stochastic process. In particular, we point out when such a reformulation can be more pertinent 
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FIG. 1: In the middle, two typical paths x{sj3K) = Xa + A^(s) are plotted as functions of the dimensionless time s. The solid 
green path, representing a typical trajectory with a high value of /3, exhibits large fluctuations, whereas the dashed red path, 
corresponding to small /3, stays in the vicinity of the initial and final point Xa- On the right, two local-time profiles L{x) are 
shown. The broad one (solid green) arises from the violently fluctuating path x{sPh), whereas the narrow one (dashed red) 
corresponds to the path with small fluctuations. On the left, we depict a generic potential V{x). 


than the conventional “sum over histories” prescription. Secondly, we wish to outline a heuristic construction of the 
local-time representation of Pi’s. More rigorous and explicit (but less intuitive) formulation of Pi’s over ensemble of 
local times will be presented in the subsequent Section. 

To provide a physically sound motivation for the local-time representation of Pi’s we follow exposition of Paulin 
et al. in Ref. m- To this end we first consider the diagonal elements of the Bloch density matrix, i.e., p{Xa,Xa, P) 
(often referred to as the Boltzmann density). Upon shifting x ^ x + Xa, and setting x = A^, r = s^H (A = ^/ph?JM 
is the thermal de Broglie wavelength), the PI ([^ can be reformulated in terms of dimensionless quantities s and ^(s) 
as 

I rx{l )=0 c pi 

p{Xa,Xa,l3) = Y / '^^{s)expi- / ds 

^ Jx{0)=0 I Jo 

Note in particular, that in contrast to 'Dx{t) the measure T’C(s) does not explicitly depend on /3, and thus /3-dependent 
parts in the PI are under better control. Such a rescaled representation is particularly useful when discussing large- 
and/or small-/3 behavior of the path integral in question. Path fluctuations in the potential are controlled by A oc 
and the factor /3 quantifies the significance of the potential V with respect to the kinetic term. 

For small /3 (i.e., high temperature), typical paths x{sf3fi) = a;a+A^(s) stay in the vicinity of the point Xa-, as depicted 
in Fig. [l] and therefore a systematic Wigner-Kirkwood expansion can be readily developed by Taylor-expanding the 
potential part of the action m- 

When /3 is large (i.e., low temperature), the trajectories x{sj3H) fluctuate heavily around the value Xa, and the 
potential term V dominates over the kinetic one. From the statistical physics point of view, the most important 
contribution to the low temperature behavior of the path integral should come from those paths that spend a 
sizable amount of time near the global minimum of the potential V (x). For this reason, it is important to be able to 
keep track of the time which a given path spends in an infinitesimal neighborhood of an arbitrary point x. 

Let us define, for each Wiener trajectory x(t) present in the Feynman path integral ^ the ensuing local time as 

L^{t) = f dr'6{X — x{t')) . (5) 

Jo 

Since the local time (t) is a functional of the random trajectory x(t') for 0 < r' < r, it represents a random 
variable. From the definition ^ we can immediately see that > 0 for all Af S M, fgdX L^(t) = t and that 
has a compact support. In addition, it can be proved mm that local-time trajectories are, with probability 
one, continuous curves which (similarly as trajectories in the underlying Wiener process) are nowhere differentiable. 
In Fig. we depict two examples of representative local-time trajectories. An extensive mathematical discussion of 
properties of the local time can be found, e.g. in Refs. [HIS]. 


+ (JV{Xa + A^) 


(4) 
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With the definition ([^ the potential part of the Euclidean action can be recast into form fgdX (/3h)V(X). A 
local-time representation of the Bloch density matrix p{xa,Xb, /3) is then given by 


p{xa,Xb,l3) = J'DL^W[L;P,Xa,Xb]5(JdXL^-l3ti\exp -^JdXL^V{X) 


( 6 ) 


where the PI “sum” is taken over all local-time trajectories L^, with x being the independent variable (not to be 
mistaken with the Wiener trajectory x{t)). The 5-function enforces the normalization constraint mentioned above. 
Basically, transition to the local-time description represents a change (or a functional substitution) of stochastic 
variables x{t) — )■ The weight factor W appearing in ®> can be formally written in the form 


'yV[L]f3,Xa,Xb] 


exp 



det 


/ 6L^i(3h) 

\ Sx{t) 



(7) 


It is a function of Xa, Xb (which are implicitly present in and /3, and a functional of the local time . Of 

course, these cavalier manipulations do not have more than a heuristic nature, and it is, indeed, a non-trivial task to 
determine W directly from 0. For this reason, we will in the following Section tackle this problem indirectly. 


IV. LOCAL-TIME REPRESENTATION OF PATH INTEGRALS: DERIVATION 


In this Section, we present a derivation of the local-time representation of the Bloch density matrix (IT). Initially, 
we limit ourselves to considering the case of the diagonal part, Xb = Xa, and arrive at the key result (Tm, which 
expresses the matrix elements in the ensuing Laplace picture with respect to /3. This intermediate outcome is shown 
to agree with the Sturm-Liouville theory. In the next step, we generalize the latter result to the off-diagonal elements 
{xb ^ Xa)- The inverse Laplace transform will then yield the sought local-time representation of PI [cf. Eq. (28)]. 


A. Field-theoretic representation 

It follows from the definition ([^ that the function p{xa,Xb, 13) satisfies the heat equation 

p{Xa,Xb,/3) = 0, 


dl3 2M dxl ^ 


( 8 ) 


with the initial condition p(xa,Xb,0+) = S{xa — Xb)- This is merely a Wick-rotated (t —)■ —iUP) analogue of the 
Schrodinger equation. The Feynman-Kac formula [Duniiss] then ensures that the PI ([^ can be calculated by solving 
corresponding parabolic differential equation ([^. 

In the Laplace picture Eq. (1^ takes the form 


^2 92 ^ ' 


p{Xa,Xb,E) = 5{Xa-Xb)^ 


(9) 


with p{xa,Xb, E) = d/Se p{Xa,Xb, l3). Eq. (|^ implies that p is nothing but the Green function of the operator 
E + H. With the benefit of hindsight, we represent the Green function p as a path integral over fluctuating fields — 
so-called functional integral [5]. This is rather standard strategy in Quantum Field Theory [iiin], and in our case it 
yields 




p{Xa,Xb,E) = 




( 10 ) 


where 


rX+ 




dx ip^x) 


IX- 

.x+ 




1 p(x) 


dx 


—,p'(x)^ + {V{x)+E)i;{x)^ 


( 11 ) 
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is the Euclidean action functional of the field-theoretic path integral. The super-index “E” in A indicates the shift 
in the potential V(x) by the amount E. Here, we have conhned our quantum-mechanical system within a finite box 
with X_ minjxajXb} and X^ ^ max{xa,Xb}. A real scalar field ip{x) satisfies Dirichlet boundary 
conditions '(/;(A_) = '0(A_|_) = 0 so as to ensure the validity of the operations being performed. 


B. Replica trick 


Since we will ultimately want to invert the Laplace transform to regain from p{xa, x^, E) the original Bloch density 
matrix p{xa,Xb, P), we cannot treat the denominator in ( |10[ ) as an irrelevant normalization constant (which is the 
usual practice in Quantum Field Theory) but we have to take care of its E-dependence. To this end we take advantage 
of the formula 


— = lim ab 
b D^O 


D-l 


( 12 ) 


which is a simple version of the so-called replica trick. [The usual replica-trick formula (cf. e.g., Ref. |36| ) can be 
obtained from (12) by integrating both sides with respect to b and subsequently dividing by a.] With the help of (12) 
we can rewrite (10) as a multidimensional functional integral 


p{Xa,Xb,E) = lim — / T>'j/)(a;)'0(a;a) 


(13) 


where the multiplet xjj = (i/ii,... ,iPd) is a E-component “replica” field in 1 -1- 0 dimensions, ip{Xa) ■ ip{xb) denotes 
the scalar product i’a{xa)4’a{xb), and we have rescaled the helds by a factor of y/2 in passing. The factor 1/E 

results from a PI generalization of the well know mean-value identity {xi pi) = (x • y)/E valid for any two vectors in 
E-dimensional statistically isotropic environments. 

As a side remark, note that we may now invert the Laplace transform, using the trivial identity 


POO 

/ dj3 5[p — c) = for c > 0, 

Jo 


to obtain the representation 


2 r'4’ix+)—0 

p{xa,Xb,P) = lim — / Dtp ^p{xa) ■ tpixb) Si 

D^o D J ^( X _)=0 


tp{xYdx — /3 I e 




(14) 


(15) 


Upon scaling tp — >• y/JJip, this agrees with the formula (2.9) in Ref. [IH]. Though the result (15) generalizes to arbitrary 
number of dimensions of the x-space, i.e., x G our further development will be illustrated (for simplicity’s sake) 
only on the one-dimensional case. 


C. Connection with radial harmonic oscillator 


In order to derive the weight factor >V (cf. Eq. Q) we have arrived at the representation (15) with E replica 
fields. This form is still not very transparent, and a further simplification step is needed to get rid of an explicit 
dependence of the measure on E. To this end, we note that action of a E-dimensional 

harmonic oscillator with the “time” variable x, “position” variable ip, mass hJ/M, and time-dependent “frequency” 
V{x) + E. Considering for a moment only diagonal matrix elements, Xb = Xa, spherical symmetry in the replica field 
space allows to reduce the path integral (13) to its radial part. Due to the boundary conditions, %p{X-) = 'i/j(A+) = 0, 
only the zero-angular-momentum (s-wave) contribution is non-vanishing (generally a weighted sum over radial Pi’s 
with different angular momenta would be required). This will be rigorously justified at the end of this Section. The 
corresponding radial PI representation for (13) reads 


p{Xa,Xa,E) = lim 


rv{x+)=v+ 


—r lim (p-'nA.y 
D—>0 D^l(D) >-0 


'v(X-)=V- 


'Diri{x) rf{xa) e _ 


(16) 


Here, the radial part rj = \/^ of the E-dimensional replica field ip is always non-negative, i.e., ri{x) > 0; the area of 
a unit sphere in E dimensions, E(E) = 2Tr^^'^IT{D 12), may be replaced by its small-E asymptotic form H(E) ~ E; 
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and ? 7 ± have been introduced to regularize the origin of the tp-space. The new action functional 






dx 


/x_ 


M {D- l){D-2,) 
h? 8?7^(x) 


(17) 


is the Euclidean action functional of the radial harmonic oscillator [ZlElllST]. It contains an additional centrifugal 
potential term (Edwards-Gulyaev or Langer term [71 [13[35]), which emerges from Bessel function Id/ 2 -i present in 
the finite sliced form of the radial PI (161. At this point we should stress that in contrast to the quantum-mechanical 
radial PI, one can use safely the asymptotic expansion for the Bessel function Id/ 2 -i (see, e.g. Ref. 


]): 




1 


-l/4)/2yj 


{\yj \ > 1, Re[yj] > 0), 


yj = {m/eh)rjrj_i , (18) 

in the Euclidean PI sliced form. Here the infinitesimal “time” slice e is related to the number of slices N via the 
relation e = hfi/N. In quantum mechanic this is a problematic step because (18) requires Re[?/j] > 0 while there 
Re[j/j] = R,e[[m/ieh)rjrj-i] = 0. 

Fortunately, the PI for radial harmonic oscillator is exactly solvable even in the case of x-dependent oscillator 
frequency. The solution reads [24] 


rVA2)=ri2 r t-X2 

{ri 2 X 2 \riiXi)D = / T>77(x) exp / dx 

Jr]{xi)-Tii I Jxi 


,2 2 M(D-1)(£I-3) 


M Gixi) 


D/2-1 


h? 


mm \ 


(F'{x 2)^ 
M G{xi) J[ 2 M\F{x2) 


exp 


-m 


G\x,) 2 

G(x,) 


The functions F{x) and G{x) are two independent solutions of the differential equation 


H + E 


y{x) = 


h?' d^ 


■E 


y{x) = 0, 


(19) 


( 20 ) 


with the initial conditions F{xi) = 0 and E'{xi) = 1, and G{x 2 ) = 0 and G'{x 2 ) = —1. In addition, the Wronskian 
W(F, G) = F(x)G'(x) — F'(x)G(x) is independent of x, as can be proved by differentiation and by using the fact that 
F and G both satisfy Eq. (20). By equating the values of IT(E, G) at points xi and X 2 , and taking into account the 
initial conditions for F and G, we find a useful identity F{x 2 ) = G{xi). 

Now, the PI in (16) can be sliced at point Xa, and expressed as 


poo 

p{Xa,Xa,E) = / dr]a{V+X+\r]aXa)D'nl{r]aXa\'n-X-)D ■ 

Jo 


( 21 ) 


The limits in Eq. ( [I^ ) are readily carried out with the help of the asymptotic formulas 
Id/ 2 -i{z) « (z/2)'°/^’^r(ZI/2), and r( 2 ;) « l/z, valid for z —>• 0+. Subsequent integration over r/a 

brings (16) to form 


p(X(j, Xfj, E) 


Fl{Xa)G2{Xa) 


2M 


Fi{Xa)G'^iXa) - F{{Xa)G2iXa) 


( 22 ) 


where Fi{x) solves Eq. (l2^ with initial conditions Fi(A_) = 0 and F{(A_) = 1, and G 2 {x) solves the same equation 
with G 2 (X_|_) = 0 and G^A+) = — 1. The denominator in (22) is the Wronskian IT(Fi,G 2 ). The full derivation is 
given in Appendix [Aj 

Although rather explicit, Eq. ( [2^ is not well suited for the Laplace transform inversion, since functions F{x) and 
G{x) contain A in a non-trivial way, which, in addition, significantly hinges on the actual form of V(x). For formal 
manipulations it is still better the employ the PI representation (16). For instance, using Eq. (14) we can easily invert 
the Laplace transform to return from E back to the /3-variable, namely 


p(a;a,Xa,/3) = lim ^ lim (?7_77+) ^ 

£>->■0 V ±^0 


rv{x+)=v+ 


lrj(X_)=V- 


V-qiyX) r] {xa)S 



y^dx — /3 I e 


(23) 
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Note that we have utilized the asymptotic form n{D) ^ D which holds for D We shall see shortly that (23) can 
be straightforwardly related to the local-time PI representation of the Boltzmann density matrix. 

Let us finally comment on the higher-angular-momentum terms which, as claimed, should not contribute to the 
expression (16). For arbitrary angular momentum ^ > 0, we employ formula (19) with a slight modification D —>■ 

> 0, this goes like {rjaXa\ri-XJ)D+ 2 i oc which, multiplied by the 

ryl. That is, only the (£ = 0)-term can give a non-vanishing contribution. 


Zl-|-2^. Now, for example, in the limit ry_ 
prefactor d/ 2 ^ implies the behavior 


D. Connection with the Sturm—Lionville problem 


Consider again Eq. ([9|) and a finite interval x G [X-,X+]. The corresponding Green function of the operator H + E 
can be easily constructed (at least formally) with the help of the Sturm-Liouville theory [53]. An immediate 
consequence of the latter is that for Xa < Xb the Green function has the form 


p{Xa,Xb,E) = - 


2MF{xa)G{xb) 
1T(F, G) 


(24) 


where the functions F{x) and G(x) satisfy Eq. (20) with the initial conditions F{X_) = 0 and F'{X_) = 1, and 
G{X+) = 0 and G'{X+) = — 1, respectively. In addition, the above Green function should be symmetric due to the 
Hermitian nature of El. 


The Sturm-Liouville theory ensures that the solution to the second-order differential equation (20) is unique, when 


specifying the values of y(a;o) and y'ixo) at some point Xq. Therefore, the functions F and G must coincide with Fi 
and G 2 of Eq. (22), and the diagonal part of (24), i.e. p(xa,Xa,F), reduces to expression (22). This is an important 


consistency check of our representation (16). 


E. Extension to off-diagonal matrix elements 


Let us now generalize the PI representation to the full Bloch density matrix, i.e., we wish to include also the off- 
diagonal matrix elements, Xb ^ Xa- If we go back to the replica representation (13), we realize that the requirement 
Xb ^ Xa spoils rotational symmetry in the replica field space, and thus precludes straightforward reduction to a 
radial path integral. Instead of refining the reduction procedure, we simply make a guess, which, as we prove in 
Appendix coincides with the well-established Sturm-Liouville Formula (241. Our guess is based on mathematical 
results presented in m- In particular, we claim that the extension of the representation ( |16| ) to off-diagonal matrix 
elements should read 


p{xa,Xb,E) = lim lim {r]_r]+)~ / Vpix) r]{xa)ri{xb) e~ , 

D^O v±^0 Jn(X.)=v. 


ri{X-)=ri- 


with the action functional 


rX. 


Aflv] ^ A^lv] + 


+ , M A(x) 
ax 


ix_ 


lA 8ri‘^(x) 


where is defined in Eq. (11), and A{x) is a piecewise constant function 

. . _ |-1 for X e [Xa,Xb] 

\ {D — 1){D — 3) otherwise. 


(25) 


(26) 


(27) 


At this point we can invert the Laplace transform with the help of Eq. (14). As a result, we obtain the sought 
local-time PI representation of the Bloch density matrix ([^, namely 


p(Xa,X6,/3) = lim ^ lim (?y_ry+) 2 
D-fO V‘‘ v±^0 


rv(.X+)=rj+ 


/rj{X-)=T]- 


Vt] r]ixa)r]ixb)d\ 


rfdx — /3 e 






(28) 


/x_ 


Here, integrations over ?y(x) run from 0 to -l-oo, i.e., the paths ?y(x) are non-negative. Comparing this result with 
the anticipated heuristic form ([^, we can identify ^^{x) = L^{j5h)/h. Representation (28) allows us to identify the 



















weight factor Q with 


yVoimP^Xa^Xb] = ^(? 7 _? 7 +)- 


. [ /•^+ 
r]{xa)r]{xb) exp < - y dx 


_,2 , M A{x) 
h? 


2M 




(29) 


Contrary to expectation, the right-hand-side of this expression does not depend on 13. Sub-index D in Wd indicates 
that the weight factor must be regularized when we pull it out of the PI (28). By analogy with quantum mechanics 


one can represent (28) in the discretized time-sliced form. In such a case the weight Wd would be a product of terms 
involving the Bessel functions /_d/ 2 -i, if A{x) = {D — 1){D — 3), or Jg, if A{x) = —I (see Ref. [H]). 

Last but not least, expressions (28)-(29) indicate that the square root of is (at least from physicist’s point of 


view) more convenient variable to describe local-time trajectories than alone. 


V. FUNCTIONALS OF THE LOCAL TIME 


Formula (28) provides a way of rewriting the PI (§ in terms of the local time. In this Section, we consider more 
general scenario, in which the initial path integral is of the form 


rx(l3h)=Xi r 2 ^ I'Ph 

F{xa,Xb,f3) = / I?a;(T)F[L] exp </ dr 

J x{ 0 )—Xa I ^ Jo 


M 2 

—X +V{x) 


(30) 


where F is an arbitrary functional of the local time L^{/3ti), which itself is (as seen in Section III) a functional of the 
paths x{t). Relation (30) represents a local-time analog of the Feynman-Matthews-Salam formula 


To bring it into more manageable form, we may observe that for any X, the action of in the PI (30) can be taken 




over by the functional derivative —hd/SV{X), acting on the exponential. This becomes transparent after rewriting 
the potential part as f^^dr V(x(t)) = J^dX V{X)L{X). The entire functional F[L] can be therefore pulled out of 
the path integral, which then allows to write 


F{Xa,Xb,(3) = F 




p{Xa,Xb,l3) . 


(31) 


When we employ the local-time representation of PI for p{xa, Xb, /3) (cf. Eq. (28)), each functional derivative 
—hS/SV(x) will produce the term hr]^(x). In such a way F can be written as 


Fi^Xfi., X 5 , / 3 ) 

2 rvix+)=7j+ 

= lim ^ lim ( 77 - 77 +) 2 / Vp'q{xa)'n{xb) 5\ I p^dx -/3 ] F[hp^]e 

D-s-O v±^0 Jr,iX.)=r,. Jx. J 







(32) 


where, strictly speaking, the functional F’[- • ■ ] is regularized in such a way that it depends on only for X € 
[A_, A+], and X± are sent to ±00 only at the end of the calculation. 

First, let us make the simple observation that the formula (32) reduces toj|28|) for the choice F[L\ = 1. One of the 


most important mean values of a local-time functional, as evaluated with Eq. (32), is the mean of exp(— dXL^j(X)) 
which gives the moment-generating functional. The local-time moment structure is particularly pertinent in various 
perturbative expansions, including low- and high-temperature expansions (see Section VI). Another important exam¬ 


ple, namely the case of a one-point distribution function will be discussed in the following subsection. 


In passing we should note, that should we have started from (15) and repeated the above procedure, an analog 
of Eq. (32) for higher-dimensional spaces, x € could be easily obtained. This would include the H-dimensional 


replica field 1/7 in the d-dimensional Euclidean configuration space. 


A. Example: One-point distribution function at the origin 

Simple, though quite important consequence of Eq. ( |32| ) is that it readily provides the 7V-point distribution functions 
of the local time. This is achieved when we set F[L] = 11^=1 ^(-^^" ~ Ln). In order to see what is involved let us 
now illustrate the calculation for N = I (with Li = L). Our discussion will be greatly simplified by considering only 
a free particle (i.e., V{x) = 0) that starts and ends at the origin, i.e., Xa = Xb = 0. This corresponds to a stochastic 
process known as Brownian bridge. Our goal is to derive the one-point distribution function, denoted p(L; ^5), of the 
























9 


local time at X = 0. We define p{L; jS) by Eq. (30) with F\L] = 5{L^ — L), and calculate it from the representation 
(32) as follows. 

In the Laplace picture, p{L\E) = df3e~^^p{L-, (3), the path integral (32) can be sliced a.t Xa = Xb = 0 so that 


2 ID 

= lim lim (??-??+)^ / dpop'^ SiHrjl - L){r]+X+\r]o 0 )d{,Vo 0\v-X-)d ■ 


(33) 


The 770 -integration can be done easily by realizing that S{ Hp^ — L) = S{po ~ \lFjK)/2\/hL. Furthermore, the limits 
in p± and D can be carried out with the help of formulas (A2) and (|A3[) from Appendix Consequently, we obtain 


p{L; E) = exp 


M2 rF[{Q) G'M 


2M \Fi{0) GsiO) 


(34) 


where, for the free-particle case, ^ 1 ( 0 :) = sinh[-\/2M£'/Ii2(x — X-)]/yj2ME/K^, and Gz{x) = su\h[-\j2 ME/K^{Xj^ — 
X )y^2ME/f9, as one can straightforwardly verify. In the limit X± —>■ ± 00 , Eq. (^34| reduces to 


p{L-,E) = e-V^^^y 


and its inverse-Laplace transform yields 


p(L;/3) = 


^exp (-gg) 

y/2TTMI3^/E^ 


(35) 


(36) 


We stress that p{L]j3) thus obtained is, in fact, the (unnormalized ) joint probability density for stochastic events 
“x( 0 ) = 0 x{l3h) = 0” and “L° = L”. By Bayes’ theorem of the probability calculus, the desired conditional 
probability density p[L° = L\x{0) = 0-^ x[(}h) = 0] is obtained from (36) by dividing p(L; (3) by the Brownian-bridge 
probability density p[x(0) = 0 ~^a;(,dl 7 ) = 0], which is (omitting again normalization) {2TTj3ti^ (see, e.g.. Ref. [4]). 
Normalization factors mutually cancel in the fraction and we arrive at 


p[L^ = L\x{Q)=Q~^x{l3K)=Q\ = 


h^L 


exp 


2/3M J 


I3M 


(37) 


which is clearly normalized to 1. One can proceed along the same lines also in more complicated higher-dimensional 
{N > 1) cases. Our result agrees with the one found through other means in Ref. [I3]. 


VI. ASYMPTOTIC BEHAVIOR OF THE BLOCH DENSITY MATRIX 


A compelling feature of the local-time representation (28) is that it naturally captures both small- and large-/3 
asymptotic regimes. This should be compared with the Feynman-Kac PI representation ([^, which is typically suitable 
only for the small-/3 (i.e., large-temperature) analysis. The latter is epitomized either by WKB approximation [6l[7j or 
Wigner-Kirkwood expansion m- In the large-/? (small-temperature) limit, the spectral representation of the Gibbs 
operator, |(/)„) {(jn], reduces the Bloch density matrix to the ground-state contribution 


p{Xa,Xb,P) 




i’LMi’gsixb) 


(38) 


that is not evident from the Feynman-Kac PI representation |27j . In connection with Eq. (38) it is useful to remind 
that in d = 1 the discrete bound states can all be chosen to be real |32] , so that the Bloch density matrix is real and 
symmetric and can be written in the form 


p{Xa,Xb,P) 


'^^'^1pn{Xa)yn{Xb) 6 ijgsiXaHgsiXb) ■ 

n—0 


(39) 


Let us first comment on the small-/? regime of the local-time representation, assuming Xb = Xa for simplicity. This 
case was discussed in detail in our previous article |21j . There one should first Taylor-expand the potential V (x) 
around the point Xa, and then expand the exponential part containing the structure f rj^(x)0(x — Xa)dx, where 

O(x-Xa) = /3J2 J -Xa)r. (40) 

m /0 
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After the term jX is factored out of the integral, the individual summands of the ensuing series are of the 

form (32) with the potential V(x) = 0, and functional F[L] oc 1^- latter can be related to the power 

expansion in j3 presented in through the equality of representations (30) and (|32|). The whole Bloch density 


matrix (containing also off-diagonal elements) can be treated similarly in a full analogy with Ref. pT] , 

The large-/3 expansion of Eq. (28) can be conveniently studied after rescaling 77 —>■ in which case we can write 


2 1 ^ rv{x^)=r,+ 

p{xa,Xb,/3) = lim lim ( 77 _? 7 +)^ / Vr]{x) r]{xa)v{xb) 

D^O ri±^0 




X si 


rfdx — 1 exp < — 


dx 


lx. 


lx. 


/3n^ 


M A{x) 


if’'"++ i 3 h^ 8,2 j 


('C-\-ioo 


dn 




i,To ' SJ{Xa)SJ{xb) 


'Dr]{x) 


X exp < —13 


dx 


lx. 


-b E {x)^^ - Kri^] +K 


nX^ 


X exp ■ 


lx. 


M A{x) 

7 - + Jv 


, 1=0 


(41) 


where c is an arbitrary real number. With the method of images [a EH Eg we can rewrite the radial PI involved as a 
superposition of two genuine one-dimensional PI ’s[S7] 


r’i(^+)=')+ 


V- 


h{X.)=r]~ 


* 77 ( 31 ) exp I-/3 ( 77 |iJ| 77 ) - k(( 77 | 77 ) - 1) 


exp 


lx 


M A(x) 




ra(x+)=rj+ 

Jr]{X.)=Tj. 

X exp 

— cos{ttD/2) 


Vr]{x)expl^-(3 {r]\H\r]) - k {{rjlTj) - 1) 
M A{x) 


lx. Ph? 877 ' 

rv{X+)=r]+ 


lri(X.) = -7]. 


X> 77 (a;)exp |-/3 (Tyj^lTy) - k (( 77 I 77 ) - 1) -(J|| 77 |)| 


X exp 




/x_ 




(42) 


where Dirac’s notation was employed. A few comments are in order about the right-hand-side of the above relation. 
First, it should be noticed the presence of the parity-even terms (J|| 77 |) in Pi’s. Second, Pi’s differ by their respective 
Dirichlet boundary conditions. Third, the super-index “i?” was used to stress the restricted nature of the fluctuations 
in the radial PI measure while the measure without “i?” represents a usual one-dimensional PI measure, i.e. 


/ of.2 yv/^fv-i „oo / 0^2 \N/‘iN-l „oo 

P*,M=„to(^) n I ( 2 ^) <“'*> 

^ ^ fc=l ^ ^ k=l 

Here “=” denotes De Witt’s “equivalence” symbol [30]. Finally, the correct time-sliced form of the exponential with 
the centrifugal potential is (cf. e.g.. Refs. 13 EH) 


exp 


exp 


f^+ , M A{x) 
lx. ^ PfP 8772 


N 


= lim TT \ 27t 
N^oo J-J- V 


rX. 


+ , M A{x) 
dx 


j3h? d>rf- 


N 


_Ph? 77fc77fc-] 

M e 


Ak exp - 


Ph^ dfcdfc-i 
M e 


Ai 


Phi^ ijkiik-i 
M e 


■ y TJ I I Pfl^ \VkVk-l\ 7 

= J™ 




M 


(44) 
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with 


Afc = A(xfc) = 


-{D-l){D-3) fOT Xk € [Xa,Xb] 
1 otherwise, 


and (see, e.g.. Refs. [371113]) 


ippix) = 


sin(7rp) 


— [I-p{x) - Ip{x)] = e^\ —Kp{x). 


(45) 


(46) 


\Ip and Kp are the modified Bessel functions of the first and the second kind, respectively.] In cases when x ^ 1 
meaning that \r]kr]k-i\ S> e (e.g., “typical situation” for very fine time slicings) the asymptotic form of ipp{—x) ^ 
1 + (1 — 4p^)/8a; + 0{l/x^) holds. With help of the preceding asymptotic behavior one obtains 


N 




D-2 \ — 


IvkVk-i 

M e 


N 


n 

k^l 
N 

n 






exp 


8|77fe77fc_i|Afe 

M Afc 


fX 


I3h? 8|77fc?7fc_i 

M A(x) 


£ + 0(e") 


(47) 


Potential singularities of the integral at 77 = 0 can be regularized, e.g., by a principal value prescription. Unfortunately, 
the formula (471 cannot be directly used in our case because the boundary values rj- and rj^ are arbitrarily close to 
zero, and hence the assumed asymptotic behavior for is not fulfilled. This situation can be rectified by factorizing 
out the problematic boundary points as 


N 




D-2 1 — 




\rikr]k-i 

M e 


■A, 


TT \2Me) ''' 




(D-l )/2 


N-1 


n 


fc=2 


M 


Ph? 8|77fc?7fc_] 


■ e + C>(£^) 


(48) 


Here we have utilized the asymptotic form tpp{—x) ~ (—a;/2)^’+^/^ r{—p)/y^ + valid for a; <C 1 and p < 0 

(p^z-). 

The passage from the radial PI (41) to the ordinary (1-dimensional) PI brings about an important advantage, 
namely, one can perform the WKB approximation. In particular, one can use the Laplace’s formula of the asymptotic 
calculus [371133] 


dtf{t)e^^p[-Pg{t)] = 


277 


Pa" (to) 


/(to) exp [-Pg{to)] + O 


exp [-Pg{to)]\ 

p3/2 




(49) 


with to being a solution of g'it) = 0 (provided g{t) has a smooth absolute minimum at the interior point t = to ( 7 ^ 
± 00 )). In case of need, the full asymptotic expansion can be found, e.g., in Ref. [33]. In our case the customary PI 
substitution 


27r 


1 —^ lim 


Pa" (to) N —^00 V 27tsM 


ph? 


det„ ( -£"VV + e" {V{x) - Eo) 


- 1 - 1/2 


(50) 


should be employed [Slzlisi]. Here Eq is the ground state energy and the difference operators (lattice derivatives) 
are defined as [7] 


Vg{x) = -[r 7 (a;-|-£) - vix)] , V? 7 (a;) = - [f]{x) - g{x - e)] 

£ £ 


(51) 


with VV = VV being the Hermitian operator on the space of “time-sliced” functions with vanishing end points, 
i.e., ^{xn) = vixo) = 0. The ensuing determinant det„(- • •) in (50) can be conveniently computed by means of the 
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Gelfand-Yaglom method [H] (or shifting method 0 )- There the determinant det„(- • •) = D{xn)/£ where D{x) is 
finite in the e —> 0 limit and fulfills the differential equation 


dx"^ 


2M 

If 


iV{x) - Eo) 


D{x) = 0, 


(52) 


with the Cauchy conditions D{X-) = 0 and D'{X-) = 1. The solution D{x) can be written as a linear combination 
of two independent solutions of (52). One is immediate, namely the WKB solution rjQ^x). A second solution can be 
constructed with the help of d’Alambert’s formula directly from r]o{x) in the form (cf. e.g., Ref. [53]) 


Voix) 


fX. 


VoiO 


(53) 


It is simple to verify that (53) is again a solution of (52). The Cauchy conditions in the Gelfand-Yaglom method 

)rm 

(54) 


uniquely fix the constants in front of the two independent solutions so that the full solution D{x) acquires the form 

, , r d^ 


D{x) = r]o{x)rio{X_) = r]o{x)ri_ 


lx. vUO 

As a result we obtain in the large N (i.e., small e) limit 


vm ■ 


det. 


-e^XX + e^^{V{x)-Eo) 


e^O r]+TJ_ 




lx. 


dx 

rilix) ■ 


(55) 


The great advantage of the Gelfand-Yaglom method is that it does not require any detailed knowledge of the spectrum 
of the operator whose determinant is being computed. It also shows that the determinant in (50) diverges linearly 
with N in the large N limit. 


We now substitute for /(to) in (49) the functional expression 


exp 


r* dx ^ 


/x_ 


fdff StjI 


(56) 


where, r]o{x) comes from the path that minimizes the functional {rj\H\ri)—K {{r]\r]) — 1). According to the Rayleigh-Ritz 
variation principle (see, e.g.. Refs. [31] [35]), such a function r]{x) is the ground-state wavefunction of the Hamiltonian 
H, i.e., r]Q{x) = ipgsix) with k = Eq. With a hindsight, we have used the latter value of k already in Eqs. (50), (52) 


and (55). Notice that the stationary point in k is real but the integration contour in k is parallel to the imaginary 


axis. Both reality and positivity of ri{x) pose no restriction in the Rayleigh-Ritz principle, because the ground state 
can always be chosen real and positive [33]. Putting it all together, we get the leading large-/3 behavior in the form 


p{Xa,Xb,P) = e 


- 


'>Pgs{Xa)lpgs{xb) , 


(57) 


as it should be (cf. Eq. @). 

We conclude the discussion of low-temperature expansion by noting that the Rayleigh-Ritz variation principle states 
that all eigenvalues and (normalized) eigenvectors of E[ come from stationary solutions of {ri\E[\r]) — k ((?7|?7) — 1), and 
conversely [55]. In the spirit of the WKB approximation one should sum over all path integrals evaluated about 
all stationary solutions. It is, however, only the ground state configuration {ipgs{x), Eq} that acquires the global 
minimum and which gives the larges contribution to the WKB approximation. This fact was implicitly used in our 
preceding reasonings. Should we have included also other stationary solutions we would recover higher order terms 
in the spectral expansion of the Bloch density matrix (39). 


VII. CONCLUSION AND OUTLOOK 

In this paper, we have derived the local-time PI representation of the Bloch density matrix. We have shown that 
the result obtained, apart from being of interest in pure mathematics (stochastic theory, Sturm-Liouville theory, etc), 
can serve as a useful alternative to the traditional Feynman-Kac PI representation of Green functions of Fokker- 
Planck equations. Furthermore, by analytically continuing the result back to the real time via the inverse Wick 
rotation, /3 —>■ it/h, one obtains the local-time PI representation of quantum-mechanical transition amplitudes, i.e.. 
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matrix elements of the evolution operator From the physicists’ point of view, perhaps the most important 

application of local-time Pi’s lies in statistical physics, and namely in the low and high temperature treatments of 
the Bloch density matrix. This is because in conventional Pi’s only a very tiny subset of paths gives a relevant 
contribution in these asymptotic regimes. In particular, the high-temperature regime of the Boltzmann density 
function p(a;, a;, /3) is dominated by paths that spend a sizable amount of time in the vicinity of the point x. Similarly, 
the low-temperature regime is controlled by paths with a large local time near the global minimum of the potential. 
Here we have exemplihed the conceptual convenience of the local-time formulation by providing a generic analysis of 
the low-temperature behavior of the Bloch density matrix. Our formulation proved to be particularly instrumental in 
obtaining the correct asymptotic behavior (known from the spectral theory) which is otherwise notoriously difficult to 
obtain within the Feynman-Kac PI framework iiE]. As a byproduct we have uncovered an interesting connection 
between a low-temperature PI expansion, the Gelfand-Yaglom formula and the Rayleigh-Ritz variational principle. 

In order to further reinforce our analysis, we formulated a local-time analog of the Feynman-Matthews-Salam for¬ 
mula which is (similarly as its QFT counterpart) expedient in number of statistical-physics contexts. The prescription 
obtained was substantiated by an explicit calculation of a one-point distribution function of the local time. 

It appears worthwhile to stress that our local-time representation (with its build-in replica field trick) is in its 
present form applicable only to one-dimensional quantum mechanical systems. With a hindsight we reflected this fact 
already in our choice of the incipient PI ([^ where we assumed t S K and a; S IR. Though one may easily proceed up 
to Eq. (15) without any restriction on the value of d (in fact, Eq. (15) is valid for any a; S with d > 1), a further 
progress in this direction is hindered by the fact that the replica fields depend on a d-dimensional argument x, and 
thus the PI in (15) can no longer be regarded as a quantum mechanical PI (i.e., PI over fluctuating paths). In effect, 
we cannot use existing mathematical techniques of the PI calculus (e.g., transformation of Pi’s to polar coordinates), 
that we have employed to get the radial PI (16). The issue of the extension of our local-time PI representation to 
higher-dimensional configuration space is currently under active investigation. 
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Appendix A: Off-diagonal matrix elements 


In this Appendix we show that the representation (25) reduces to the well established result (24) of the Sturm- 
Liouville theory. Since p(xa, Xf,, E) is symmetric in Xa and Xb, we will assume, without loss of generality, that Xq < Xf,. 
The path integral in (25) can be expressed via Eq. (19) as 


dijadrib {r]+X+\ribXb)Dr]bir]bXb\r]aXa)2r]a{r]aXa\r]-X-)D ■ 


The limits in Eq. (25) can be carried with the help of the asymptotic formulas Id/ 2 -i{z) 
and r( 2 :) « 1/z, valid for z —>■ 0. We obtain 


1 — E> 

{'naXa\'n-X_)D -!■ 


/ -2 'n/^pxnf td F[(x^) 2 

( \ tixp ^ 2M Fi{xc,)'la 


\2MF^{Xa] 




D-^-O 




exp - 


F[{Xa) 

2MF,iXa) 


Va 


(Al) 

(z/2)^/2-Vr(D/2), 


(A2) 


where E’i(x) satisfies Eq. (20) with initial conditions F'i(A_) = 0 and F{(A_) = 1, and similarly, we find 


1 i-p 

—??+" {r]+X+\r]bXb)D 


? 7 +— >^0 1 


G'^jxb) 2 
2M Gsixb)'^^ 


(A3) 
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where Gs^x) satisfies Eq. (20) with initial conditions G 3 (X_|_) = 0 and G 3 (X+) = — 1. Formula (25) then reduces to 


2h‘^ f 

p{x,,x,,E) = 


f VaPb 

G2{Xa)^\M G2{Xa) 


, , VaVb . 

dr]adr]b ^ ^ h 


_h^ f W{G2,F,)pl W{G^,F2Wb \ 

^ 2M \Fi{Xa)G2{Xa) F2{Xb)G3{Xb) ) _ 


(A4) 


where F 2 {x) and G 2 {x) satisfy Eq. ( [^ with initial conditions F 2 {xa) = 0 and F^^Xa) = 1, and G 2 {xb) = 0 and 
G' 2 {xb) = —1 , respectively. The Wronskian W{F,G) = F{x)G'{x) — F'{x)G{x) is independent of x, as discussed in 
Section IV C[ and antisymmetric, i.e., IF(F, G) = —W{G,F). 

Wronskians kF(G 2 ,F’i) and IF(G 3 ,F’ 2 ) assume a particularly simple form when evaluated at points Xb and Xa, 
respectively, due to the initial conditions satisfied by G 2 and F 2 . We find kF(G 2 ,F’i) = Fi{xb) and kF(G 3 ,F’ 2 ) = 
Gz{xa)- Rescaling r]a —>■ y/G 2 {Xa)M/h?r]a, r]b —>■ y^F 2 {xb)M/h?r]b, and using the relation F 2 {xb) = G 2 {xa) we obtain 


_ 2M f°° 

p{xa,Xb,F) = -^G2{xa) J dpadrjb VaVb loiVaVb) exp 
The integrations are readily performed using the formula |34j 


F\{xb^ 2 G^i^Xa) 
2F^)^‘^~ 2G^) 


Vb 




yielding 


p{Xa,Xb,F) = 


2M Fi{Xa)Gz{Xb)G2{Xa) 

Fi{xb)Gi,{Xa) - Fi{Xa)G^{xb)' 


To prove equality with (24), we only have to show that 

Ei(a;o)G3(xt,) - Fi(xh)G3{xo) = G2(a:a)IF(Ei, G3). 


(A5) 


(A6) 


(A7) 


(AS) 


This is done by realizing that G 2 (x), being a solution of the second-order linear differential equation (20), can be 
uniquely composed as a linear combination of two other solutions Fi{x) and G 3 (a;), 


G2{x) = 


Fi{x)G^{xb) - Fi{xb)G3{x) 
IE(Fi,G3) 


(A9) 


Indeed, thus defined G 2 satisfies the initial conditions G 2 {xb) = 0 and G 2 {xb) = —1. 
We conclude that 


p{Xa,Xb,F) = - 


2M Fi{Xa)G3{Xb) 

W(Fi,G3) ’ 


(AIO) 


which coincides with the Sturm-Liouville result (24). 
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